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Abstract. In this paper, The central black hole masses 
and the Doppler factors are derived for PKS 0528+134, 
PKS 0537-441, 3C279, PKS 1406-074, PKS 1622-297, 
Q1633+382, Mkn 501, and BL Lacertae. The masses ob- 
tained are in the range of (1 - 7) xl0 7 M Q and compared 
with that obtained with the Klein-Nishina cross section 
considered (Dermer & Gehrels 1995). If we considered only 
the Thomson cross section, the masses are in the range of 
2.6xl0 6 M Q - 2xl0 11 M Q . The masses obtained from our 
method are less sensitive to the flux than those obtained 
from Dermer & Gehrels (1995) method. The masses ob- 
tained from two flares (1991 and 1996 flares) of 3C279 
are almost the same. For 3C279 and BL Lacertae, view- 
ing angle, 9, and Lorentz factor, T, are estimated from the 
derived Doppler factor and the measured superluminal ve- 
locity. For 3C279, 9 = 10°. 9 - 15°.6, V = 2.4-14.4 for 6 = 
3.37; 9 = 8°.45 - 9°. 7, T = 2.95-11.20 for S = 4.89; For 
BL Lacertae, 9 = 25° - 29°.4, T = 2.0-4.0. 
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1. Introduction 

One of the most important results of the CGRO/EGRET 
instrument in the held of extragalactic astronomy is the 
discovery that blazars (i.e., flat-spectrum radio quasars- 
(FSRQs) and BL Lac objects) emit most of their bolomet- 
ric luminosity in the high 7-rays ( E > 100 MeV ) energy 
range. Many of the 7-ray emitters are also superluminal 
radio sources (von Montigny et al. 1995). The common 
properties of these EGRET-detected AGNs are the fol- 
lowing: The 7-ray flux is dominant over the flux in lower 
energy bands; The 7-ray luminosity above 100 MeV ranges 
from less than 3 x 10 44 erg s _1 to more than 10 49 erg s _1 
(assuming isotropic emission); Many of the sources are 



strongly variable in the 7-ray band on timescales from 
days to months (Mukherjee et al. 1997), but large flux 
variability on short timescales of < 1 day is also detected 
( see below ). Some correlations between 7-ray and the 
lower energetic bands are discussed (see Dondi & Ghis- 
ellini 1995; Fan 1997; Fan et al. 1998a; Miicke et al. 1997; 
Xie et al. 1997; Zhou et al. 1997). These suggest that the 
7-ray emission is likely from the jet. 

Various of models for 7-ray emission have been pro- 
posed: Namely, (1) the inverse Compton process on the 
external photons (ECS), in which the soft photons are 
directly from a nearby accretion disk ( Dermer et al. 
1992; Coppi et al. 1993 ) or from disk radiation repro- 
cessed in some region of AGNs ( e.g. broad emission 
line region) (Sikora et al. 1994; Blandford & Levinson 
1995); (2) the synchrotron self-Compton model (SSC), 
in which the soft photons originate as synchrotron emis- 
sion in the jet ( Maraschi et al. 1992; Bloom & Marscher 
1992, 1993; Zdziarski & Krolik 1993; Bloom & Marscher 
1996; Marscher & Travis 1996); (3) synchrotron emission 
from ultrarelativistic electrons and positrons produced in 
a proton-induced cascade (PIC) (Mannheim & Bicrmann 
1992; Mannheim 1993; Cheng & Ding 1994). TeV radi- 
ations are observed from 3 X-ray-selected BL Lacertae 
objects (XBLs): Mkn 421 ( Punch et al. 1992); Mkn 501 
(Quinn et al. 1996), and IES 2344+514 (Catanese et al. 
1998). But there is no consensus yet on the dominant 
emission process (see 3C273 for instance, von Montigny et 
al. 1997); for PKS 0528+134, the lower and higher states 
can be fitted by different models (eg. Bottcher & Collmar 
1998). 

The 7-rays are produced at a distance of ~ 100i? s 
(Hartman et al. 1996), 205R g (Xie et al. 1998), and hun- 
dreds of Schwarzshild radii (Ghisellini & Madau, 1996; 
Celotti & Ghisellini 1998). We think that this distance 
is an important parameter, which can be used to con- 
strain the mass of the central black hole. In this paper, we 
will use it to derive the central black hole mass and the 
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Dopplcr factor for some blazars with short timescales. The 
paper is arranged as follows: In section 2, we estimate the 
mass of the central black hole and the Doppler factor; in 
section 3, we give some discussions and a brief summary. 
Ho = 75 km s _1 Mpc -1 , and qo = 0.5 are adopted 
through out the paper. 



2. The mass of the central black hole and the 
Doppler factor 

2.1. Data 

Some 7-ray loud blazars have been observed several 
times with EGRET while three XBLs have been ob- 
served to show TeV radiations. In the following section, we 
present the 7-ray loud blazars with available 7-ray vari- 
ation timescales. Because the variability timescale corre- 
sponds to different variation amplitude for different source 
and/or different observation period, we use the doubling 
timescale, ATd = (Tj n j t , a j/AT)AT, as the variability 
timescale. 



2.1.1. PKS 0528+134 

PKS 0528+134, z = 2.07 (Hunter et al. 1993), is one of 
the most luminous examples of blazars. It is observed by 
EGRET, COMPTE and OSSE aboard the CGRO (see 
Hunter ct al. 1993, McNaron-Brown et al. 1995, Mukherjee 
et al. 1996, Collmar et al. 1997; Sambruna et al. 1997). 

During the period of 16-30 May 1991, the source 
showed F(>100 MeV) = ( 1.0 ± 0.2) xl0~ 6 photon cm~ 2 
s _1 with photon spectral index a 7 = 2.56 ± 0.09. 

During 23-29 March 1993, F(>100 MeV) = (0.23+0.12 
- 3.08 ± 0.35) xl0~ 6 photon cm" 2 s _1 ; with a photon 
spectral index a 7 = 2.21+0.10. In the 1993 observation, 
a variation of order 100% over a timescale of <~ 2 days 
was detected (see Wagner et al. 1997), which suggest a 
doubling time scale of ATd = lday. 

During August, 1994, F(>100 MeV) = (0.32 ± 0.1) 
xlO -6 photon cm~ 2 s _1 ; with a photon spectral index a 7 
= 2.70. 

There is a clear evidence that the spectrum becomes 
harder when the 7-ray flux increases. 



2.1.2. PKS 0537-441 

PKS 0537-441, z = 0.896, a candidate of gravitational 
lens (Surpi et al. 1996 ), is a violently variable object ( 
Fan & Lin 1998). The 7-ray flux varies from (1.83+0.91) 
to (8.98+1.45) xl0~ 6 photon cm" 2 s" 1 (Mukherjee et al. 
1997). A flare of a factor of -3 from 0.35 to 2.0 xl0~ 6 
photon cm _2 s _1 over a time scale of <~ 2 days can be seen 
from Fig. 3 in Hartman's paper (Hartman 1996). AT/3 = 
16 hrs. 



2.1.3. 1253-055, 3C279 

3C279 is a well known member of OVV subclass of blazars. 
it is perhaps the prototypical supcrluminal radio source 
(Moffct ct al. 1972); and the first quasar detected at the 
energies of >1 GcV with EGRET/CGRO. The simultane- 
ous variability in X-rays and 7-rays (> 100 MeV) suggest 
for the first time that they are approximately cospatial ( 
M c Hardy 1996). The 7-ray flux varies from 1.28 to 28.7 
xl0~ 6 photon cm~ 2 s _1 (Mukherjee et al. 1997). Two 7- 
ray flares were detected (see Kniffen et al. 1993; Hartman 
ct al. 1996; M c Hardy 1996; Wehrle et al. 1998). 

The 16-28 June 1991 flare showed: F(>100 MeV) = ( 
2.8 ± 0.4) x 10~ 6 photon cm~ 2 s _1 with a photon spectral 
index a 7 = 1.89+0.06. A variation of a factor of 4 over 2 
days was obtained. 

The January- February 1996 flare showed ( see 
McHardy 1996; Wehrle et al. 1998), F(>100 MeV) = ( 
11.0 ± 1.) xl0~ 6 photon cm~ 2 s _1 with a photon spec- 
tral index a 7 = 1.97+0.07. During this flare, a variation of 
a factor of 4^5 in a day is observed, AT = 6 hrs (Wehrle 
et al. 1998). 

No obvious spectral index variation has been detected 
when the flux varied. 



2.1.4. PKS 1406-074 

PKS 1406-074 has been detected to vary with the 7-ray 
flux being in the range of 1.54 to 12.76xl0~ 6 photon cm -2 
s _1 (Mukherjee et al. 1997). During its 7-ray flare, a flux of 
F(>100 MeV) = ( 5.5 ± 1.4) x 10~ 6 photon cm" 2 s" 1 with 
a photon spectral index a 7 = 2.04+0.15 and a doubling 
timescale of shorter than 16 hours has been obtained (see 
Wagner et al. 1995). 



2.1.5. PKS 1622-297 

For PKS 1622-297, 2=0.815, we have very little informa- 
tion in lower energy bands. But it is one of the most 
luminous objects in the 7-ray region. A peak flux of 
(17+3) xl0~ 6 photon cm~ 2 s" 1 (E>100 MeV) and a flux 
increase by a factor of 2 in 9.7 hours were observed ( Mat- 
tox et al. 1997). 



2.1.6. Q1633+382, 4C 38.41 

Quasar 1633+382, z = 1.814, is an LPQ (P opt = 2.6%, 
Moore & Stockman 1984). During 1992 November 17 - De- 
cember 1 period, it was detected to show a flux of F(>100 
MeV) = ( 0.30 ± 0.06) xl0~ 6 photon cm" 2 s" 1 with a 
photon spectral index a 7 = 1.87+0.07. The flux varied a 
factor of 1.5 within 24 hrs, ATp = 16 hrs, while the spec- 
tral index did not change. The 7-ray luminosity is at least 
two orders of magnitude larger than the maximum ever 
observed in any other band (see Mattox et al. 1993). 
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2.1.7. B2 1652+399, Mkn 501 

Mkn 501, 2 = 0.033, together with other two XBLs are 
three known TeV 7-ray sources detected by the Whip- 
ple group ( see Quinn et al. 1996; Catanese et al. 1997a; 
Samuelson et al. 1998; Kataoka et al. 1998). 

During 1995 March-July period, Mkn 501 was observed 
to show a flux of F(>300 GeV) = ( 8.1 ± 1.4) xl0~ 12 
photon cm~ 2 s _1 with a photon spectral index a 7 = 2.2. A 
variation of a factor of 4 over one day is also detected, ATp 
= 6 hrs. The upper limit corresponds to a flux of F(>100 
MeV) = 1.5 xl0~ 7 photon cm~ 2 s _1 (sec Quinn et al. 
1996). During the 1996 multiwavclcngth campaign, Mkn 
501 was detected with EGRET a flux of F(>100 MeV) = 
(0.32 ± .13) x 10~ 6 photon cm~ 2 s _1 with a photon index 
of 1.6±0.5 (see Kataoka et al. 1998). During 1997 April 9- 
19 observation, Catanese ct al. (1997a) obtained F(>300 
GcV) = ( 40.5 ± 9.6) x lO" 11 photon cm" 2 s _1 , a 7 = 2.5, 
the April 9-15 flux corresponds to a flux of F(>100 MeV) 
< 3.6 xl0~ 7 photon cm~ 2 s _1 . 

The TeV observations show that the spectrum softens 
when the source brightens. 

2.1.8. 2200+420, BL Lacertae 

2200+420 is the prototype of BL Lacertae class. It is 
variable in all wavelengths (see Fan et al. 1998b, 1998c; 
Bloom ct al. 1997; Bottchcr & Bloom 1998; Madejski et 
al. 1998). A 14- year period was found in the optical light 
curve (Fan et al. 1998b). During 1995 January 24 - Febru- 
ary 14, BL Lacertae showed a flux of F(>100 MeV) = ( 
40 ±12) xl0~ 8 photon cm~ 2 s _1 with a photon spectral 
index a~ f = 2.2+0.3, the up limit flux in higher energy is 
F(>300 GeV) < 0.53 xlO" 11 photon cm" 2 s" 1 ( Catanese 
et al. 1997b). During 1997 January 15/22 observation pe- 
riod, it was detected a flux of F(>100 MeV) = ( 171 ±42) 
xlO -8 photon cm~ 2 s _1 with a photon spectral index a 7 
= 1.68+0.16 and a dramatic factor of 2.5 increase within 
a timescale of 8hrs, ATq = 3.2 hrs. Besides, simultaneous 
optical and 7-ray flares were observed ruling out external 
scattering models (sec Bloom et al. 1997). 

The observations from the object show that the spec- 
trum of BL Lacertae hardens when the 7-ray flux in- 
creases. 



2.2. The central black hole mass and the Doppler factor 

The objects discussed here show variability time scale of 
hours to days. The variability could be directly related to 
shock processes in a jet, far from the accretion disk ( we 
thank Dr. S. D. Bloom to point out this for us). If we take 
the variability timescale as the measurements of the size, 
R, of the emission region, then the R in the jet obeys to 
the inequality, 



R < cAT D 



where c is the speed of light, S the Doppler factor, 2 the 
redshift of the source, and ATq, in units of second, the 
doubling time scale. 

For an object with a mass M, the Eddington limit 
gives (Frank, King & Raine 1985) 



00 M 
L Edd . w 1.26 x 10 38 (— -)ergs- 

Mq 



(2) 



So, we have that the intrinsic luminosity, L vn " of a source 
with a mass of M should satisfy L m < LEdd.- 

In the relativistic beaming frame, the observed lumi- 
nosity is L ob - = (5( 4+a ) L m - , a is the energy spectral index, 
which follows that 



L ob. < d (4+ a ) L 



Edd. 



(3) 



Ghiscllini & Madau (1996) obtained that the 7-rays 
are emitted within the BLR region, which is 10 17-18 cm 
far from the central source. Hartman et al (1996) obtained 
that the 7-rays are produced at a distance of <~ 100i? 9 . Re- 
cently, Celotti & Ghiscllini (1998) argued that the 7-rays 
are from a region of some hundreds of Schwarzschild radii 
from the center. From our previous paper, a distance of 
205R g is obtained for the 7-rays from Mkn 421. In the 
sense of the theory of accretion (Sunyaev 1975). When 
R < 200R g , the electrons in the accretion flow become ul- 
trarelativistic. On the other hand, the mixture of relativis- 
tic electrons and nonrelativistic protons has an adiabatic 
index 7 < |, with such an adiabatic index the transition 
to supersonic accretion regime is possible in the region 
R < 200R g (Sunyaev, 1975). So, the 200i? s is perhaps an 
important critical point. If we assume that the 7-rays are 
from this place then relations (1), (2), and (3) give 



M „ <5 

= 5 x 10 2 

M Q 1 + 2 



AT, 



D 



A(5+«) 

L ob - < 6.3 x 10 40 - -AToergs- 1 

(1 + z) 



(4) 



(5) 



It is a common property of the EGRET-detected 
AGNs to show that their 7-ray flux is dominant over 
the flux in lower energy bands but this is not always 
the case (Mukherjee et al. 1997). For PKS 0528+134 
and 3C279, their 7-ray luminosity, L 1 , is 0.80ifc o ;. and 
0.5Lboi., respectively (see Sambruna et al. 1997; Hartman 
et al. 1996). Because we consider the flare states of the se- 
lected objects, we can take the 7-ray luminosity to stand 
for half of the bolometrical luminosity approximately, i.e. 
L 7 ~ 0.5 x Ltd.- So, we have 

x(5+a) 

L 7 < 3.15 x 10 40 - -ATnergs- 1 

(1 + z) 

which gives 



(6) 



(1 + *) 



cm 



(1) 5>[ L " {l + z) ](*&> 

v ; - L 3.15 x 10 40 ergs- 1 AT D i 



(7) 
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So, from the available L 7 and ATd, we can obtain the 
central black hole mass and the Doppler factor from rela- 
tions (4) and (7) and a = a 7 — 1, a 7 is the photon spectral 
index. They are shown in table 1, in which, Col. 1 gives the 
name; Col. 2, the redshift; Col. 3. the flux F(>100MeV) 
in units of 10~ 6 photon cm~ 2 s , a is the uncertainty; 
Col. 4, the photon spectral index, cu 7 = 2.0 is adopted 
for 0537-441 (see Fan et al. 1998a); Col. 5, the doubling 
time scale in units of hours; Col. 6, 7-ray luminosity (as- 
suming isotropic emission) in units of 10 48 erg s _1 ; Col. 7, 
Doppler factor estimated from equation (7); Col. 8, the 
central black hole mass in units of 1O 7 .M ; Col. 9, the 
central black hole mass estimated from Dermer & Gchrcls 
(1996), thereafter D&G, in units of 10 7 M Q ; Col. 10, the 
mass estimated directly from Eddington limit in units of 
1O 1O M 

3. Discussion 

3.1. Mass 

From the high 7-ray luminosity (assuming isotropic emis- 
sion) and Eddington-limit, one can derive the central black 
hole mass expression, 



M w > 



Ln 



1.26 x 10 48 ergs- 



(8) 



where, L<r is the bolometric luminosity for emission in 
the Thomson region, M10 is the central black hole mass 
in units of 1O 1O M0. The derived masses are as high as 
10 n M o for some 7-ray loud blazars, PKS 0528+134, PKS 
1406-074, and PKS1622-297 for instance (see Col. 10 in 
table 1). However, for high energy 7-ray emission, Klein- 
Nishina effects must be considered. D&G considered the 
effect and obtained an expression for the black hole mass, 
i.e. their equation (16b), 



Mi CN > 



3ird 2 L (m e c 



F(ei,e v 



2 x 1.26 x lO^ergs- 1 1 



-Zn[2e«(l + z)] (9) 



where F(si,s u ) is the integrated photon flux in units of 
10~ 6 photon cm~ 2 s _1 between photon energies si and 
e u in units of 0.51 lMeV. For the objects considered here, 
Mf- N is obtained and shown in Col. 9 in table 1. Table 
1 shows that the masses obtained from our consideration 
and those estimated from D&G method are acceptably 
similar except for 1622-297 and two low redshift BL Lac 
objects (Mkn 501 and BL Lacertae). For 1622-297 our 
value is about 7 times less than that estimated from D&G 
method. If we adopt the flux density 2.45 x 10 -6 photon 
cm -2 s _1 for 1622-297 instead of the peak value as did 
Muhkerjee et al. (1997) and Fan et al. (1998a), then the 
isotropic luminosity is 3.87 x 10 48 erg s _1 . This luminosity 
suggests that the Doppler factor and mass obtained from 
relations (4) and (7) are respectively 5.01 and 2.4IM7, and 
the mass estimated from D&G is then 3.6IM7. The both 
masses are quite similar in this case. For Mkn 501 and BL 



Lacertae, our value is much greater than that estimated 
from D&G method. For 3C279, our results show that the 
estimated central black hole masses, 4.74M7 and 3.43 M7 
for 1991 and 1996 flares respectively, are almost the same, 
while the masses estimated from D&G method are I.92M7 
and 7.53M7 for 1991 and 1996 flares respectively. Table 1 
(also see relation (9)) shows that the mass obtained from 
D&G method is sensitive to the flux, variable flux gives 
different mass for a source. The mass obtained from our 
method does not depend on the flux so sensitively For 
3C279 1991 and 1996 flares, masses obtained from our 
consideration are almost the same while those obtained 
from D&G method show a difference of more or less a 
factor of 4. But our method depends on the timescales (see 
relation 4). Since we only considered the objects showing 
short timescales ( hours) in the present paper, the masses 
obtained are in a range of (1 ~ 7) xlO 7 M . 

To fit 3C279 multiwavelength energy spectrum corre- 
sponding to 1991 7-ray flare, Hartman employed an ac- 
creting black hole of 1O 8 A# , our result of 4.74 x 1O 7 M 
is similar to theirs. 

3.2. Beaming factors 

To explain the extremely high and violently variable lu- 
minosity of AGNs, beaming model has been proposed. 
In this model, the Lorentz factor, T, and the viewing 
angle, 8, are not measurable, but they can be obtained 
through the measurement of superluminal velocity, (3 app , , 
and the determination of Doppler factor, 5, which are 
related with the two unmcasurablc parameters, T and 
6, in the forms: [3 app = n-£!t"% > r = 7=^ > and 



5 = (r(l — PinCosQ)) 1 . So, r and 8 can be obtained from 
the following relations: 



r = 



0iw + s2 + 1 



26 



= tan ( 



— 1 { ^Papp 

~W 2 



+ 6 2 - 1 



From our previous work (Fan et al. 1996), we can get 
superluminal velocities for 3C279 and BL Lacertae. When 
the superluminal velocities and the derived Doppler factor 
are substituted to the above two relations, we found that: 
For 3C279, T = 2.4-14.4 and 8 = 10°. 9 - 15°. 6 for 6 = 
3.37; and T = 2.95-11.20 and 8 = 8°.45 - 9°. 7 for 6 = 
4.89. For BL Lacertae T = 2.- 4. and 8 = 25° - 29°.4. 

To let the optical depth (r 77 ) be less than unity, 
Doppler factor in 7-ray region has been obtained for some 
objects by other authors. (5 > 7.6 for Q1633+382 ( Mattox 
et al., 1993); <5 > 6.3 - 8.5 for 3C279 1996 flare ( Wehrlc 
et al. 1998) and 5 > 3.9 for 3C279 1991 flare ( Mattox et 
al. 1993), 5 ~ 5 is also obtained by Henri et al. (1993); 
5 > 6.6 ~ 8.1 for PKS 1622-297 (Mattox et al. 1997). Our 
results in table 1 are consistent with those results. 
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Table 1. Mass and Doppler factor for 7-ray loud blazars 



Name 


z 


F(a) 


Ct-y 


AT D 


-^48 


8 


M 7 


M 7 KN 


Ml 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0528+134 


2.07 


3.08(0.35) 


2.21 


24. 


18.4 


4.96 


6.97 


22.98 


14.6 


0537-441 


0.894 


2.0(0.4) 


2.0 


16. 


3.01 


3.83 


5.82 


3.48 


2.39 


1253-055 


0.537 


2.8(0.4) 


2.02 


12. 


1.34 


3.37 


4.74 


1.92 


1.06 


1253-055 


0.538 


11.(1.) 


1.97 


6. 


5.75 


4.89 


3.43 


7.53 


4.56 


1406-074 


1.494 


5.5(1.6) 


2.04 


16.21 


23.7 


5.57 


6.52 


23.6 


18.81 


1622-297 


0.815 


17.(3.) 


1.87 


4.85 


26.9 


6.97 


3.35 


25. 


21.35 


1633+382 


1.814 


0.96(0.08) 


1.86 


16. 


9.72 


5.16 


5.28 


5.74 


7.71 


1652+399 


0.033 


<0.36 


2.5 


6. 


0.0003 


0.89 


0.94 


0.001 


0.00026 


2200+420 


0.07 


1.71(0.42) 


1.68 


3.2 


0.019 


2.04 


1.09 


0.02 


0.015 



3.3. Summary 

In this paper, the central mass and Doppler factor are 
obtained for 8 7-ray loud blazars with available short 7- 
ray timescales. The mass obtained from relation (4) in 
the present paper is compared with that obtained from 
D&G method, the masses obtained from two methods are 
similar for 5 out of 8 objects. Our method is not sensitive 
to the flux as much as D&G method in estimating the 
central black hole mass. The masses obtained here are in 
a range of (1 ~ 7) xl0 7 M Q , which is from the fact that 
the time scales considered here are in a range of 3.2 to 
24 hours. For 3C279, the masses obtained from the two 
flares are almost the same. For 3C279 and BL Lacertae, 
the Lorcntz factor and viewing angle are estimated. 
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